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Abstract 

A group of transformations changing the phases of the elements of 
the single-particle density matrix, but leaving unchanged the predic- 
tions for identical particles concerning the momentum distributions, 
momentum correlations etc., is identified. Its implications for the de- 
terminations of the interaction regions from studies of Bose-Einstein 
correlations are discussed. 

PACS 25.75. Gz, 13.65.+i 
Bose-Einstein correlations. 



1 Introduction 

Bose-Einstein correlations are helpful when trying to derive properties of 
the interaction regions from the measured momentum distributions. In this 
report we point out an ambiguity inherent in such derivations pp. There are 
many ways from the data to the inferred properties of the interaction regions 
(cf. e.g. j2] and references given there). One can use density matrices, 
Wigner functions, emission functions, distances between the pairs of points 
where the identical particles are produced etc. The ambiguity seems to be 
common to all of them. 

We will not discuss here the evolution of the density matrix during the 
freeze-out period. Formally, one can introduce the assumption that all the 
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hadrons have been produced instantaneously and simultaneously at some 
time which we may choose as t=0. Then what one measures is the density 
matrix (in the interaction representation) at freeze-out. 

As is well known, the diagonal elements of the fc-particle density matrix 
in the momentum representation give, and are unambiguously given by, the 
k particle momentum distribution. On the other hand, in most models these 
elements can be expressed as symmetrized products of the single particle 
density matrix elements [Sj: 

k 

p(Pi,---,Pfc;Pi,---,Pfc) = ^Ylpii.Pj'^Ppj)^ (1) 

p j=i 

where the summation is over all the k\ permutations of the momenta p^^, . . . , p^i. 
Thus, all the momentum distributions are unambiguously determined when 
the single particle density matrix pi(pi; P2) is known. 

Our main observation ^ is that the converse is not true. Given the 
momentum distributions for all the sets of k = 1,2, .. . particles, it is not 
possible to find unambiguously the matrix pi. This trivial observation will 
be seen to have very non trivial consequences. Since the matrix pi is further 
used the derive conclusions about the interaction region, the ambiguity affects 
our capacity for making unambiguous statements about such regions. 

2 Invariance group 

Consider the transformation 

Pi(p; P') - Pi.(p; P') = e^"^pVi(p; pOe-^^^^'^ (2) 

where a is any real-valued function of momentum. According to formula 
(|H) for = 1, pi is a single particle density matrix. Therefore, it must 
be hermitian and have trace one. Also pia, as seen from its definition, is 
hermitian and has trace one. Consequently, it can be substituted for pi on 
the right-hand side of formula (0). This introduces on the right-hand side, for 
every pi a term e'"^P»^ and a term e"*"^^'-' which cancels it. Thus, the diagonal 
matrix element on the left-hand side does not change. Experimentally, the 
substitution of pia for pi is invisible. The transformations from pi to pia 
form a local (in momentum space) U(l) invariance group. 

There are several quantities related to the single particle density matrix in 
the momentum representation and yielding information about the interaction 
region. In order to get the space distribution of the sources one can use 
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the diagonal elements of the single particle density matrix in the coordinate 
representation 



pi(x;x) = jd'K^^ e^'^'Xpsp'), (3) 

where 

K=^(p + p'); q = p-p'. (4) 

In order to get an approximate phase-space distribution one can use the 
Wigner function 

W,{K,X) = /(|^e^"V(p;p'), (5) 
or the emission function 

Pi(p;p') = Jd^XSiK,X)e^'^''. (6) 

In the last two formulae X = ^{x + x') and the approximation consists in 
interpreting K and X as the energy-momentum and space-time position of 
the particle. This approximation is not always good, but in a well-defined 
sense |3] this is the best one can have without contradicting the principles 
of quantum mechanics. In the last formula the four-vector K has the same 
value on both sides of the equality. Since the momenta p and p' are on 
mass shell, for any q the value of go is fixed by the condition Kq = 0. 
Thus, it is not possible to invert the Fourier transform and to express the 
emission function S in terms of the density matrix pi. In fact, there is 
an infinity of different emission functions, corresponding to various off mass 
shell continuations of a given on mass shell function pi. We conclude that 
the ambiguities when using the emission function formalism are more severe 
than when using Wigner functions. 

The point is that the transition from pi to pia, which has no effect on 
the momentum distributions, changes the functions pi, Wi and S and con- 
sequently changes the conclusions concerning the interaction region. In the 
following section we will illustrate this fact by some examples. 
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3 Examples 



The class of transformations Q is very rich. We will just discuss three 
simple examples. Consider 

a(p) = b ■ p =^ C({p) — tt(p') = b ■ q, (7) 
where b is any vector. This gives 

Pi(x;x) = JdKj^^ e^''-(-+^)pi(p;p'). (8) 

The interaction region gets shifted by b. Similarly, replacing b and p by 
four-vectors and using the emission function one can generate an arbitrary 
shift in space-time. This result is of little interest. It is obvious that the 
momentum distributions do not depend on where and when the experiment 
was done. 

As our second example consider 

"(p) = ^cp^ =^ a(p) - a(p') = cK ■ q, (9) 
where c is any real number. Then 

This time the shift is proportional to K with a proportionality coefficient 
which is unconstrained by momentum measurements. The corresponding 
distribution in space is 

pi,(x;x) = jd'KW,^{K,^) (11) 

Suppose now that for c = there are no position-momentum correlations. 
Then for each K the interaction region occupies the same portion of space 
and when averaged over K coincides with that for any given K. For c 7^ 0, 
the interaction regions corresponding to different values of K are shifted 
with respect to each other and the averaged size of the interaction region 
gets bigger. Its increase with respect to the situation at c = depends 
on the value of |c|. Since this is unconstrained by the data on momentum 
distributions, the radius can be made as large as one wishes. 
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For instance, for the Gaussian 

pi(p;p'' 



2\3 



exp 



2A2 



one obtains Piq(x;x) also Gaussian with root mean square width 



:i2) 



(13) 



In order to see that this broadening is due to the averaging over K it is 
enough to have a look at the corresponding Wigner function 



exp 



2A2 



[X + cK) 
2i?2 



(14) 



(27ri?A)3 

The conclusion from this example is that without further assumptions one 
can at best obtain a lower limit for the radius of the interaction region. In 
practice, everybody uses, more or less consciously, a model which supplies 
the necessary additional assumptions. The problem how to choose the best 
model among those which give exactly the same fit to all the data is an 
interesting open problem. In order to show that this problem is not purely 
academic let us consider the following example from the literature. 

Some models assume correlations between position in space-time 

and energy-momentum of the type^ 



A convenient notation is 



A 



(15) 



(16) 



Note that here Mt is a temporal-longitudinal variable. Let us assume that 



S = SuSt 



(17) 



where 



5- 



T 



exp 



24 



[Kt - XXrf 
25| 



This is a possible quantum-mechanical rendering of the transverse part of 
the classical relation (fT3j) . The only information about S\\ important for our 

"'^ Sometimes called Hubble flow 
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purpose is that it depends neither on K.^ nor on Xj^. The distribution of 
sources in the transverse plane is obtained by integrating St over K^. The 
result is a Gaussian with constant root mean square width tt- On the other 
hand St can be rewritten in the form 



^7 



exp 



where 



^1 



2RI 



fx. 



T 



2m 



(19) 



Rd 



Tt 



Rd — fJ'R<t>', 



t6t 



Mt 



t>T = Tt 



This can be considered as a transform with a(p) = and 

TTfl 



5T(l + /i2) 



1^ Kt 
1 + /i2 5r ■ 
(20) 



(21) 



of 



5*, 



aT 



exp 



26^ 



2Ri 



(22) 



Performing the integration of So/r over Kt one finds for the distribution of 
sources in the transverse plane a Gaussian with root mean square width 



which exhibits the familiar decrease of the transverse radius with the trans- 
verse mass as reported by so many experimental papers. Let us summarize 
the situation: if our prejudice is in favor of a Hubble fiow as interpreted in 
0, we conclude that the transverse radius does not depend on M^; if our 
prejudice is against correlations between position in the transverse plane and 
transverse momentum, we conclude that the transverse radius decreases with 
increasing My; experimental data on momentum distributions will not help 
us to decide which of these two prejudices is the right one. 

Our last example is one dimensional. It could e.g. apply to one transverse 
component. We choose the Gaussian density matrix (fT^ and 



a{p) 



3a3 



P 



a{p) - a{p') 



3a3 



q K' + 



q- 



(24) 



^This is an approximation, because c depends on Mt, it is , however, good enough for 
our quahtative discussion; compare W, where the fuU calculation can be found 
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A simple calculation |T] yields 

n 

Sn 



exp 



2A2 



where Ai{. . .) is the Airy function and 



B 



At{A), 



(25) 



A = aX + 



ijj 



B 



10 



A 



12 



X = X 



(26) 

For large values of a the emission function is almost Gaussian. A numerical 
calculation shows that a = 2/Ris already large enough. For smaller values of 
a, however, at negative X the emission function develops big wiggles, oscillat- 
ing between positive and negative values. Its shape in the positive X region 
also significantly changes This example shows how the transformations 
discussed in the present report can lead to changes of the interaction region 
which are much more complicated than just momentum dependent shifts. 



4 Conclusions 

The experimental data about momentum distributions tell us little about 
the interaction regions, unless additional assumptions are made. The usual 
recommendations: reduce the experimental errors, include more particle cor- 
relations etc. are not enough. Exactly the same fit can be obtained from 
widely different models, differing in these additional assumptions and giving 
conflicting information about the interaction region. The caveat for model 
users is that only some of the assumptions of the model are being tested by 
comparison with the data, while others, which may be very important for 
drawing inferences about the interaction region, are unconstrained by the 
data. 
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